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We use the (W-formalism to investigate the non-Gaussianity of the primordial curvature per- 
turbation in the curvaton scenario for the origin of structure. We numerically calculate the full 
probability distribution function allowing for the non-instantaneous decay of the curvaton and com- 
pare this with analytic results derived in the sudden-decay approximation. We also present results 
for the leading-order contribution to the primordial bispectrum and trispectrum. In the sudden- 
decay approximation we derive a fully non-linear expression relating the primordial perturbation to 
the initial curvaton perturbation. As an example of how non-Gaussianity provides additional con- 
straints on model parameters, we show how the primordial bispectrum on CMB scales can be used 
to constrain variance on much smaller scales in the curvaton field. Our analytical and numerical 
results allow for multiple tests of primordial non-Gaussianity, and thus they can offer consistency 
tests of the curvaton scenario. 
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I. INTRODUCTION 

Most inflationary models give rise to a nearly Gaus- 
sian distribution of the primordial curvature perturba- 
tion, (. Deviations from an exactly Gaussian distribu- 
tion are conventionally given in terms of a non-linearity 
parameter, /nl PI- The prediction for the non-linearity 
parameter from single-field models of inflation is related 
to the tilt of the power spectrum, /nl ~ n — 1 01 which is 
constrained by observations to be much less than unity. 
In principle, measurement of /nl would give a valuable 
test of the inflation, but unfortunately such a tiny non- 
Gaussianity is likely to remain unobservable. The cur- 
rent upper bound from the WMAP three-year data [|[ is 
|/nl| < 114 while Planck is expected to bring this down 
to I /nl I ;$ 5 [ij which is still orders of magnitude larger 
than the prediction for single-field inflation. 

On the other hand, multi-field models of inflation 
can lead to an observable non-Gaussianity. One well- 
motivated example is the curvaton model 0, [H, Q : I n ad- 
dition to the inflaton </> there would be another, weakly 
coupled, light scalar field (e.g., MSSM flat direction 
@j 3) j curvaton \i which was completely subdominant 
during inflation. The potential could be as simple as 
Q V = ^M 2 cj) 2 + im 2 x 2 , where the energy density of 
<f> drives inflation. At Hubble exit during inflation both 
fields acquire some classical perturbations that freeze in. 
However, the observed cosmic microwave (CMB) and 
large-scale structure (LSS) perturbations can result from 
the curvaton instead of the inflaton. If the inflaton mass 
M is much less than 10 13 GeV then perturbations due to 
the inflaton are much smaller than 10~ 5 . 
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After the end of inflation the inflaton decays into rela- 
tivistic particles ( "radiation" ) , the curvaton energy den- 
sity still being subdominant. At this stage the curva- 
ton carries an isocurvature (entropy) perturbation. The 
entropy perturbation between radiation and curvaton is 
given by S x = 3(C X — Cr)- Observations rule out purely 
isocurvature primordial perturbations [l(| [Hj], but, so 
long as the curvaton decays into radiation before pri- 
mordial nucleosynthesis, the entropy perturbation can be 
converted to an adiabatic one. This also requires that 
all the species are in thermal equilibrium and that the 
baryon asymmetry is generated after the curvaton de- 
cays. If any of these conditions are not met then the 
curvaton could also leave a residual isocurvature per- 
turbation [l^ . but in practice the amplitude of isocur- 
vature modes are severely constrained by current data 
0, EH EH an d f° r simplicity we will not consider this 
possibility in this paper. 

As the Hubble rate, if, decreases with time after in- 
flation, eventually H 2 < m 2 , and the curvaton starts 
to oscillate about the minimum of its potential. Then 
it behaves like pressurelcss dust (with density inversely 
proportional to volume, p x cx a~ 3 ) so that its relative en- 
ergy density grows with respect to radiation (p r cx a~ 4 ). 
Finally, the curvaton decays into ultra-relativistic parti- 
cles leading to the standard radiation dominated adia- 
batic primordial perturbations [46| . However, this cur- 
vaton mechanism may, from the initially Gaussian cur- 
vaton field perturbation, 8\*, create a strongly non- 
Gaussian primordial curvature perturbation, £. The non- 
Gaussianity is large if the energy density of the curvaton 
is sub-dominant when curvaton decays. Since the ampli- 
tude of the resulting perturbation depends on the model 
parameters (such as the curvaton mass m and decay rate 
r), the observational bounds on non-Gaussianity provide 
important constraints on model parameters. 
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The objective of this paper is to calculate the proba- 
bility density function (pdf) of the primordial curvature 
perturbation in the curvaton model. Since, in the early 
universe, all today's observable scales are super-Hubble 
scales after inflation, we take advantage of the separate 
universe assumption [15 , throughout the calculations 
and employ the so-called 5N- formalism [l8|, [2(| ■ This 
allows us to determine the pdf fully non-linearly (not just 
up to second or third order in the initial field pertur- 
bations) so that it will carry all the information about 
non-Gaussianity. 

Generally we can expand any field 



(fi = ip + 



n=l 



We take the background field to be spatially homoge- 
neous, <p(t), and we will further assume that the first- 
order perturbation, 6iip(t,x) is a Gaussian random field, 
consistent with what we expect from the linear evolu- 
tion of initial vacuum fluctuations. Thus the higher-order 
perturbations, 5 n ip for n > 1, will describe non-Gaussian 
perturbations of any field. 

The primordial perturbation can be described in terms 
of the non-linear curvature perturbation on uniform- 
density hypersurfaces [2~J ] 



C(t, x) = SN(t, x) + 



p(t,x) 



p(t) 



dp 



P 



(2) 



where 5N is the perturbed expansion, p the local den- 
sity and P the local pressure. We expand the curvature 
perturbation as 



C(t,x) =Ci(t,x) + V-c„(t,x) 



(3) 



where the pdf of £i is Gaussian as it is directly propor- 
tional to the initial Gaussian field perturbation, but the 
higher-order terms give rise to a non-Gaussian pdf of the 
full £. The non-linearity parameters /nl and ^nl are 
defined by 



C = Ci + ^/nlCi 2 + ^SnlCi 3 + ■ 



or, equivalently, 



C2 — ^/nlCi 



c 3 



54 
25 



5nlCi 3 



(4) 

(5) 
(6) 



The numerical factors 6/5 and 54/25 arise because in 
linear theory the primordial curvature perturbation £ is 
related to the Bardeen potential on large scales (in the 
matter-dominated era, md), <E>// m d = (3/5)£i, which im- 
plies [Emu! 



We are specifically interested in non-linear quantities 
and, as it is C n ot <f>u that is non-linearly conserved for 
adiabatic perturbation on large scales [111, 0, [H, [2(| [13] , 
we will take Eqs. ([5]) and ^ as our fundamental defini- 
tion of the primordial parameters /nl and (/nl, respec- 
tively. 

If we write the primordial power spectrum as 

(C(ki)C(k a )) - (2^) 3 P(fc 1 )<5 3 (k 1 + k 2 ) , (8) 

then the leading order contributions to the bispectrum 
and (connected part of the) trispcctrum arc given by 

(C(k x )C(k 2 )C(k 3 )) 

= (2 7 r) 3 J B(k 1 ,k 2 )^(k 1 +k 2 +k 3 ), (9) 

(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) 

= (2 7 r) 3 T(k 1 , k 2 , k 3 )5 3 (k 1 + k 2 + k 3 + k 4 ),(10) 

where 

B(k!,k 2 ) = (6/5)/ NL [F(fc 1 )P(fc 2 ) + 2pcrms] j (ll) 
r(k 1; k 2 ,k 3 ) = 

(18/25)/ NL [F(fc 1 )P(fc 2 )P(|k 1 +k 2 |) + 23 perms] 
+ (54/25). g NL[ J P(fci)P(fc 2 )P(fc 3 ) + 3pcrms] . (12) 



(7) 



The first term appearing in Eq. (|12p . which gives the 
dependence of the trispectrum on second-order pertur- 
bations, and hence /nl, was given in [2fij . But there is 
also a term dependent upon the third-order perturbation, 
and hence <?nl, which appears at the same order, and has 
a different dependence upon the four wavevectors. 

Previous estimates of non-Gaussianity in the curvaton 
scenario have been based on expansions up to second or- 
der in the curvature perturbation. We will go beyond pre- 
vious analyses and calculate the contribution of the third 
order term in Eq. (j4|) to the trispectrum. We compare our 
analytic expressions for /nl and <?nl in the sudden- decay 
approximation [20L |29| with numerical results where we 
include the gradual decay of the curvaton, transferring 
energy from the curvaton to the radiation. Indeed using 
our numerical code we are able for the first time to give 
the full probability distribution for the primordial curva- 
ture perturbation in both the sudden-decay and the non- 
instantaneous decay case. We will calculate the skewncss 
(third moment of the pdf) and kurtosis (fourth moment 
of the pdf) as well as higher moments of the fully non- 
linear probability distribution function. 

This paper is organised as follows. In Sec. [TT] we re- 
late the curvaton curvature perturbation £ x to the ini- 
tial field perturbation 8\X at the beginning of the cur- 
vaton oscillation. Then, in Sec. IIII| we derive in the 
sudden-decay approximation a non-linear equation that 
relates the primordial curvature perturbation £ (when 
curvaton has decayed) to the curvaton curvature pertur- 
bation at the beginning of curvaton oscillation, £ x , (or 
to the Gaussian curvaton field perturbation 5\* at hori- 
zon exit). We write down the full solution of this equa- 
tion in the Appendix A. In Sec. Mil we continue solv- 
ing this equation order by order, and deriving the non- 
linearity parameters /nl and (?nl in the sudden-decay 
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approximation. In Sec. IIVI we describe our fully non- 
linear numerical approach. We compare the numerical 
non-instantaneous decay results for /nl and <?nl to the 
sudden-decay approximation. In Sec. [V] we calculate the 
pdf of £ both in the sudden-decay approximation and in 
the non-instantaneous decay case. We show the skewness 
and kurtosis of the pdf as a function of curvaton model 
parameters. In Sec. I VII we add one possible complica- 
tion to the analysis. Namely, the variance of curvaton 
field perturbations ((Sx) 2 /x 2 ) ori smaller than observ- 
able scales is not directly constrained. However, the large 
variance leads to large non-Gaussianity (see e.g. (30l.[3~lj) 
so that the observational bounds on non-Gaussianity set 
an upper limit to this small-scale variance. We derive a 
quantitative equation that relates /nl and the variance, 
and use this equation with WMAP third year bounds. 
Finally, in Sec. IVIll we present concluding remarks. 



where we used the relation 5\\ = g' ' 5\X* and wrote g 
S (X.) and <?(«) = d n g( X )/d X n \ M . 

Substituting (TJHD and JT7D into © wc obtain 



E 



9 Six 
9' X 



Order by order, we have from (fT5|) 
Sip x = m 2 g5ix, 



S 2 p x 
and from HI 



9£l 

n l2 



2!o9" , 99'" 



(Six) 3 



(19) 

(20) 
(21) 

(22) 



II. NON-LINEARITY OF THE CURVATON 
PERTURBATION 

When the curvaton starts to oscillate about the mini- 
mum of its potential, but before it decays, the non-linear 
curvature perturbation on uniform-curvaton density hy- 
persurfaces is given by [2l[ 

Px(t) 3 PX 

Hence, the curvaton density on spatially-flat hypersur- 
faccs is 



C x (t,x)=5JV(t,x)+ / 



Px\SN=0 



(14) 



Assuming the curvaton potential is described by a 
quadratic potential about its minimum, the energy den- 
sity is given in terms of the amplitude of the curvaton 
field oscillations 



Px = 



(15) 



We expect the quantum fluctuations in a weakly cou- 
pled field such as the curvaton at Hubble exit during in- 
flation, Sx*, to be well described by a Gaussian random 
field (see e.g. [H, HH [34| ) . Hence we will write 



X* 



Six* 



(16) 



with no higher-order, non-Gaussian terms. 

Non-linear evolution on large scales is possible if the 
curvaton potential deviates from a purely quadratic po- 
tential away from its minimum [35ll3q|. Thus, in general, 
the initial amplitude of curvaton oscillations, x, is some 
function of the field value at the Hubble exit; x = <?(x*)- 
(The curvaton potential is in any case virtually quadratic 
sufficiently close to the minimum.) Thus we have during 
the curvaton oscillation 



P\ 
P\ 



1 2-2 

5 my 



|m 2 



<Lin=l nix 



( 9 §ix\ 

\g' x J 



(17) 
(18) 



Cxi 

C X 3 



2S_ix 

3 X ' 

2 ^ g" 

2\ 2g' 2 



1 9 2 g'" 



c 3 



(23) 
(24) 
(25) 



Here and in what follows, we omit the bar from g and 
simply denote it by g. 

Using (|23j) and ([24]) . we can express the second-order 
skewness in terms of the effective non-linearity parameter 
for the curvaton perturbation, analogous to Eq. ([S]), 



x 

NL 



1 



99 

g' 2 



(26) 



Hence we find /? 



, NL — —5/4 for the curvaton £ x in the 
absence of any non- linear evolution (g" — 0). If the cur- 
vaton comes to dominate the total energy density in the 
universe before it decays, so that £ = ( x , then this is the 
generic prediction for the primordial /nl in the curvaton 
model, as emphasised by [201 ] . 

From the third-order term (|25[) we obtain a contri- 
bution to the trispectrum of the curvaton perturba- 
tion, described by a non-linearity parameter analogous 
to Eq. ©, 



■9NL - 10 I 1 o „/2 + o 



25 



12 



3gg" , lg 2 g"' 



2g<2 2 g 



/:', 



(27) 



III. SUDDEN-DECAY APPROXIMATION 

Most analytic expressions for the primordial density 
perturbation in the curvaton scenario assume the instan- 
taneous decay of the curvaton particles. In this section 
we will derive an equation for the non-linear curvature 
perturbation, and then use it to find the non-linearity 
parameters /nl an d 5nl m this sudden- decay approxi- 
mation. 
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In the absence of interactions, fluids with a barotropic 
equation of state, such as radiation (P r = p r /3) or the 
non-relativistic curvaton (P x = 0), have a conserved cur- 
vature perturbation plj 



dpi 



Pi + Pi(pi) 



(28) 



We assume that the curvaton decays on a uniform- 
total density hypersurface corresponding to H = T, i.e., 
when the local Hubble rate equals the decay rate for the 
curvaton (assumed constant). Thus on this hypersurface 
we have 



Pr(tdc 



Px(t 



dec 7 ^0 



P(tdc 



(29) 



where we use a bar to denote the homogeneous, unper- 
turbed quantity. Note that from Eq. @ we have £ = SN 
on the decay surface, and we can interpret C as the per- 
turbed expansion, or ll 6N" . Assuming all the curvaton 
decay products are relativistic, we have that £ is con- 
served after the curvaton decay since the total pressure 
is simply P = p/3. 

By contrast the local curvaton and radiation densities 
on this decay surface may be inhomogeneous and we have 
from Eq. (gSJ) 



Cr 

Cx 



C + Jln 



C+ 3 ln 



or, equivalently, 



Pr 

Px 



p r e 
Px 



4(Cr-C) 

,3(C x -0 



(30) 
(31) 



(32) 
(33) 



Requiring that the total density is uniform on the de- 
cay surface, Eq. (|29|). then gives the relation 



(i-n 



X.dec, 



Mr-0 



Mx-0 



1 , 



(34) 



where f2 x .dcc = Px/(Pr + Px) is the dimensionlcss den- 
sity parameter for the curvaton at the decay time. This 
simple equation is one of the main results of this paper. 
It gives a fully non-linear relation between the primor- 
dial curvature perturbation, which remains constant 
on large scales in the radiation-dominated era after the 
curvaton decays, and the curvaton perturbation, £ x de- 
scribed in Sec. HH in the sudden-decay approximation. 
In the limiting case where ri Xi dcc - * 1 (i-e., the energy 
density of the curvaton comes to dominate before it de- 
cays) we have ( — > ( x , but in general Eq. (|3~4"|) gives a 
non-linear relation between £ and £ x . 

For simplicity we will restrict the following analysis 
to the simplest curvaton scenario in which the curvature 
perturbation in the radiation fluid before the curvaton 



decays is negligible, i.e., ( r = 0. After the curvaton de- 
cays the universe is dominated by radiation, with equa- 
tion of state P = p/3, and hence the curvature pertur- 
bation, £, is non-linearly conserved on large scales. With 
Cr = Eq. (J34J) reads 

e 4C - [ft x ,dc C e 3C *] e c + [ft x , dcc - 1] = , (35) 

which is a fourth degree equation for X = e^. In the 
Appendix A we give the solution of this equation. Since 
we already know e 3 ^ x as a function of the initial field 
perturbation 5x*, we have now found a full non-linear 
mapping of the Gaussian perturbation S\* to the pri- 
mordial (non-Gaussian) curvature perturbation We 
can Taylor expand the solution (|A2|) to find first, second, 
and third order expressions or we can (re-)solvc Eq. (|35|) 
order by order as we do in the following subsections. 



A. First order 



At first order Eq. (|3"4"1) gives 

4(1 — f2 x ,dcc)Ci = 3f2 x ,dcc(Cxi — Ci) 
and hence we can write 

Ci = r( x i , 

where [|| 

3^ x ,dec ^Px 



4-0 



X,dec 



3p x + Ap r 



(36) 
(37) 

(38) 



*dec 



B. Second order 



At second order Eq. (|3~4"|) gives 

4(1 - J\, d ec)C2 - 16(1 - fi x ,dec)Cl 

= 3fi Xj dcc(Cx2 — C2) + 9f2 Xi dcc(Cxl — Cl) 2 

and hence, using Eqs. (|2~4")) . (f3"T|) and 

~ 3 

Jr \ g 

This gives the non-lineari ty p arameter (|5|) in the sudden- 
decay approximation [20l . l29l ] 

£ 

-.12 



<2 



99 

,/2 



cr 



5 / 99 \ 5 5r 

/NL = 47l 1+ 7^J~3~y 



(39) 

(40) 
ien- 

(41) 



In the limit r — ► 1, when the curvaton dominates the 
total energy density before it decays, we recover the non- 
linearity parameter (|26p of the curvaton 



NL 



99 
g' 2 



(42) 



On the other hand we may get a large non-Gaussianity 
(|/nl| 3> 1) in the limit r — ► [13], where we have 



/nl 



5 

ir 



1 



9j£_ 



(43) 
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Third order 



At third order we obtain from Eq. 



9 
4r 2 



9' 3 



V 



9J_ 

a' 2 



2 1 ~' 2 



+ lOr + 3r^ 



(44) 



The non-linearity parameter <?nl from Eq. ((6]) will thus 
be 25/54 times the expression in the square brackets. (As 
a consistency check we note that in the limit r — > 1 this 
result agrees with ((27)) .) 

If there is non-linear evolution of the curvaton field, 
X, between Hubble-exit and the start of the curvaton 
oscillation, such that gg" j g' 2 ~ —I, then from (f4Tj) we 
see that /nl can be small even when r — ► 0, see also 
[H, HU- However, in this case <?nl will be very large un- 
less in (j44|) the g 2 g"' /g' 3 term also cancels the 3gg"/g' 2 
term. Indeed, assuming that the g'" term is small, we find 
C 3 -» -[27/(4r 2 )]d 3 , i.e., g NL -V-25/(8r 2 ), when r -» 
if gg" j g' 2 — — 1. In this situation £3 would be of the same 
order as £1 ~ 10~ 5 , if r < 10~ 5 . Hence, even if the non- 
linear evolution of the x ne ld was such that the leading 
order non-Gaussianity, /nli was cancelled, the higher or- 
der terms could still lead to large non-Gaussianity which 
could be ruled out by observations. 

In the absence of any non-linear evolution of x field 
between the Hubble exit and the start of curvaton os- 
cillation (as in the case of truly quadratic potential) we 
would have g" = g'" = 0, so the third order result would 
be simply 



C3 

.9nl 



- + i + 10r + 3r 2 
r 2 



25 
54 



9 1 

- + - + 10r + 3r 2 

r 2 



(45) 
(46) 



It should be noted that now there is no 1 /r 2 term in this 
<7nl ■ Thus it is only at most of the same order as /nl • 
Indeed, in the limit r — > we have /nl — * 5/(4r) and 
g NL -> -25/(6r), i.e., <? nl //nl -10/3. As Ci - 10~ 5 
this means that, in this case, the third-order term (3 is 
about 5 orders of magnitude smaller than the second- 
order term £2 even when r — > 0. 



IV. NUMERICAL CALCULATION 

Although the sudden-decay approximation gives a 
good intuitive derivation of both the linear curvature per- 
turbation and the non-linearity parameters arising from 
second- and third-order effects, it is only approximate 
since it assumes the curvaton is not interacting with the 
radiation, and hence £ x remains constant on large scales, 
right up until curvaton decays. In practice the curvaton 
energy density is continually decaying once the curvaton 



begins oscillating until finally (when T > H) its density 
becomes negligible, and during this decay process ( x does 
not remain constant [13, [H[ . 

Another problem with results derived from the sudden- 
decay approximation is that the final amplitude of the 
primordial curvature perturbation, and its non-linearity, 
are given in terms of the density of the curvaton at the de- 
cay time which is not simply related to the initial curva- 
ton density, especially as the precise decay time, H ~ r, 
is ambiguous. 

In fact a more careful treatment of the continuous de- 
cay of the curvaton [3^, [38| shows that the transfer coef- 
ficient at first order, r in Eq. (|37|) . is a function solely of 
the parameter 



X M r 



(47) 



where the right hand side is to be evaluated when the 
curvaton begins to oscillate, long before it decays, and 
hence can be written as 



P 



8ir X 



2 

osc 



3Af 2 ! 



m \ !/ 2 

T 7 



(48) 



In Refs. [371 . 1 3 81 ] the resulting primordial curvature per- 
turbation in the radiation-dominated era after the curva- 
ton has completely decayed was calculated using linear 
cosmological perturbations on large scales, to give 



Ci = r(p)C x i 



(49) 



where an analytic approximation to the numerical results 
gives HU 



r(p) 



1 



1 



0.924 
1.24 



P 



-1.24 



(50) 



We find that this not only gives a good approximation 
to the amplitude of linear perturbations, but as we will 
show it can also be used to give a surprisingly accurate 
estimate for the non-linearity parameter /nl- 

In principle one could use the second-order perturbed 
field equations on large scales [3!| to evaluate C2 as a 
function of £ X 2 and hence /nl- Indeed this has recently 
been done in Ref. [40| . However a simple short-cut to 
the same result is provided by the 5 N- formalism [20I ]. 
The advantage of the JA^-formalism is that it gives im- 
mediately the results to any order one wants. Thus it 
is not necessary to repeat the calculation with more and 
more complicated perturbed field equations, if one wants 
results at higher order in perturbations. Indeed, once 
calculated, SN encodes all orders in perturbations, i.e., 
it gives the fully non-linear £. 



A. Practical implementation 

We use the evolution equations for a homoge- 
neous Fricdmann-Robcrtson- Walker (FRW) universe to 
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describe the fully non-linear evolution in the long- 
wavelength limit adopting the separate universes ap- 
proach [3, [TtJ ■ The resulting primordial curvature per- 
turbation, £, corresponds to the perturbation in the local 
integrated expansion, SN, on a final uniform-density hy- 
persurface in the radiation-dominated universe after the 
curvaton has completely decayed. We use the fully non- 
linear equations for the evolution of the homogeneous 
curvaton and radiation densities, including the gradual 
decay of the curvaton into radiation. 

Hence, our set of equations is the Friedmann equation 
and the continuity equations for curvaton and radiation 
densities. These can be written in the form I38l 



dH b 



3 + a, 



dN 
d£l x 

~dN 

da 
In 



— j — m x H mv + n r (n r — i) , 



(51) 
(52) 
(53) 



which is particularly suitable for numerical calculation. 
Here H mv = 1/H, Ct x = p x /p to t, and fl r = p r /p to t with 

Ptot = P X + Pr- 

Since the end-result does not depend on a particular 
choice of T and m, as long as we integrate far enough so 
that the curvaton has completely decayed at the end of 
calculation, we fix these to the values m = 10~ 5 Afpi and 
r = 10~ 7 m. The initial value of H lm is 1/m, since we 
start the calculation at the beginning of the curvaton os- 
cillation. After specifying the value of p, we calculate the 
initial values of Q x and Cl r . They are Q X i = (I\ffi n v) 1/ ' 2 P, 
and Q r j = 1 — f2 X j. The initial value for our integration 
variable N can be set to zero because in the absence of 
any initial perturbation in the radiation (£ r = 0) the 
initial surface is both spatially flat and has uniform en- 
ergy density and Hubble rate H = Hi = m (recall that 
from the Friedmann equation p to t oc H 2 ). We are inter- 
ested in the integrated expansion between this initial un- 
perturbed hypersurface and some final uniform-density 
surface H = Hf <C T. Then the (local) integrated ex- 
pansion between these surfaces will be just the final value 
of N = N f . 

We find that f2 x is practically zero when N > 11. From 
this we deduce that a suitable ending condition (curva- 
ton has completely decayed) is H = Hf = r/500. We 
use our modified version of an adaptive step size ode in- 
tegrator [4l[ and the accuracy parameter eps = 10~ 21 . 
We start integration with a sufficiently small step size as 
demanded for our required accuracy. Finally, when N 
starts to be of the order 11 the step trial would lead to 
-ffhw > 1/Hf. As soon as this happens we divide the 
step trial by two. We repeat this procedure until H[ nv 
obeys (1 - lCT 20 )/#/ < 4, < 1/H f . Then we save p 
and the final Nf. To find N(p) we repeat this process 
for 50000 logarithmically spaced values of p in the range 
[10~ 8 , 10 4 ]. 



B. Comparison of sudden-decay approximation 
with numerical results 



Previous studies of non-Gaussianity in the curvaton 
model have been based on the sudden-decay approxima- 
tion. In [42| we extended the calculation of the non- 
linearity parameter /nl to the non-instantaneous decay 
case and found that the sudden-decay approximation is 
indeed very accurate. Recently, a similar numerical com- 
parison was made in [40| using second order perturbation 
theory. Our results obtained using SN formalism agree 
with those of Rcf. |40|. In this subsection we describe 
our calculation of /nl [421 ] in more detail and for the first 
time perform similar studies for <?nl- 

Expanding SN we have 

C = N'6 X * + \N"(5 X *Y + ^N'"(S X ,) 3 + • • . (54) 

Comparing this with Eq. ^ we can read off — 
Qn N /Q x n^ Substituting this into © gives 



/nl 



5 N" 
6N^ 



and substituting into ([6]) we find 



25 N'' 



.9nl 



54 N 



/3 



(55) 



(56) 



As we will specify the initial conditions for our numeri- 
cal solutions by giving the value of p, defined in Eq. (|47p . 
the differentiations of N with respect to \* need to be 
converted into differentiations with respect to p. From 
we have 



d n g' d 
OX* 9 op 



Using this we find 
N' 



N = 



N = 



9 9'dN 
g dp 

, g^\ 3 d 3 N 
9 . 



2p 



dp 3 



Up 2 



9 -) 2 - 


g" 


dN 


9 J 


-— 

g 


dp 


~9'g" 




9 2 





-2p 



g_Z + 3 9Y 



dN 
dp 



(57) 



(58) 



d 2 N 
dp 2 

(60) 



Recalling the definition (|49|) of the curvature perturba- 
tion transfer efficiency at linear order, r, we find, from 
Eqs. ([23 | and ([58 ]) . 



3p 



dN 
dp 



(61) 



Once we have numerically calculated N as a function 
of p in the non-instantaneous decay case, we can calculate 
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/nl(j") and 5nl(^) by substituting (|58) 
and (|56[) . For example, for /nl we find 



ED into US 



NL 



5 

4r 



f 



gg 
g' 2 



5 d 2 N/dp 2 

6 (dN/dpf 



(62) 



Comparing to the sudden-decay result (|41j) we see that 
the first term is exactly the same. Numerical calculation 
of the derivatives of N with respect to p shows that the 
second term approaches a constant value —2.27 as r — > 
0. In the sudden-decay case it approaches a constant 
—5/3 = —1.67, see f4T]) - Thus in both cases 



-(1 

4r \ g 



gg 

12 



(63) 



when r — > 0. In the opposite limit, r — > 1, both results 
give /nl — * —5/4. So any difference between the sudden- 
decay approximation and the non-instantaneous decay 
calculation appears only at intermediate values of r, i.e., 
when the radiation energy density from the inflaton decay 
products and the curvaton energy density are of the same 
order when the curvaton decays, p r ,& C c ~ PxA^c- 

The second term in (1621) can also be written in terms 
of r. Using (JHTJ> and (JSTJ) the result ([62j) reads 



/] 



NL 



5 

4r 



I 



gg 

n l2 



5 (g/g'y 

4 r 2 



2r 



where we have used 



3P 2 



,d 2 N 
dp 2 



(64) 



(65) 



Comparing (|M|) to (|^T|) we find that in the sudden- 
decay case 



g , 



2r SD - gr SD 



3 SD ' 



(66) 



whereas in the non-instantaneous case r' must be deter- 
mined numerically. Thus one way to characterise the ac- 
curacy of the sudden-decay approximation is to calculate 
r' numerically, employing (|6f |) and (|65p . and compare it 
to the above expression (|66p in the sudden-decay approx- 
imation. 

As mentioned in the beginning of this section the rela- 
tion between r and p in the sudden-decay approximation 
is non-trivial. In the non-instantaneous decay case r(p) is 
easy to find numerically from (|6ip . In the sudden-decay 
case we can only determine rgo (p) from (|38[) if we know 

^X,dcc(p)- 

Fortunately, a short-cut to the same result is provided 
by the differential equation (|66|) . Using ((57)) we find from 



dp 



3dr SD 



3r SD 



z 'sd 



1 SD 



and hence 



poc r SD (r SD +3) 1/4 (l-r SD ) 3/4 



(67) 



(68) 
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■ Sudden decay 



0.01 
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Figure 1: To achieve the same curvature perturbation transfer 
efficiency r one needs to start from slightly different initial 
value of p in the sudden-decay case (red dashed line) than in 
the non-instantaneous decay case (black solid line). 



The constant of proportionality is not uniquely deter- 
mined by the sudden-decay approximation, and corre- 
sponds to the arbitrariness in the definition of the decay 
time H dcc ~ T. 

What we can do is to use the limiting form of the 
analytic approximation to the numerical solution (|50[) 
for small p to provide an overall normalisation for the 
sudden-decay approximation. This yields 



3I/4 



P 



rsD 



0.924(r SD + 3) 1 /4(i_ rsD )3/4 



(69) 



The above equation thus determines the value of p that 
corresponds to a given value of the linear transfer func- 
tion, tsd, in the sudden-decay approximation, and hence 



X,dec 



from Eq. 



In Fig. [T] we show p as a function 



of tsd in the sudden-decay approximation and compare 
this with the numerical non-instantaneous decay result 
for p(r). 

Our form for tsd(p) is quite different from that 
adopted by Malik and Lyth in their recent work [io| . 
They used a much simpler, but less accurate, estimate 
for r(p) in the sudden-decay approximation: 



''ML 



V 



1 



(70) 



Although tml 1 as p — > oo it fails to reproduce the 
correct linear coefficient as for p — > 0. The apparent error 
in the sudden-decay approximation for /nl reported by 
Malik and Lyth [40[ (see for instance Figure 9 in that 
paper) is primarily due to this inaccuracy in the ?"ml ijp) ■ 
In what follows we have chosen to do all our compar- 
isons of the sudden-decay approximation and numerical 
non-instantaneous decay results at common values of r. 
In other words, we have presented our results as a func- 
tion of r instead of p. After taking into account this fun- 
damental difference in thinking, the results of [1(3] agree 
with ours. 

Now we are ready for the final comparison of /nlsdMj 
derived in the sudden-decay approximation with, /nl(?") 
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Figure 2: The non-linearity parameter /nl as a function of 
curvature perturbation transfer efficiency r = Ci/Cxi- The 
analytical approximative, i.e., sudden-decay result (red dashed 
line) crosses zero at r = 0.58 and is negative for r > 0.58. 
The exact numerical result (black solid line) is negative for 
r > 0.53. Here we assume that g(x*) is linear. 



1 1 1 1 1 I 1 1 1 1 I ' 
Numerical non-instantaneous 




Figure 3: The non-linearity parameter #nl as a function of 
curvature perturbation transfer efficiency r = Ci/Cxi- The 
analytical approximative, i.e., sudden-decay result (red dashed 
line) crosses zero at r = 0.83 and is positive for r > 0.83. 
The exact numerical result (black solid line) is positive for 
r > 0.79. Here we assume that g(x*) is linear. 



calculated numerically allowing for non-instantaneous 
decay. Fig. rj shows that if / NL > 60 (r < 0.02) or 
/nl < —1.16 (r > 0.95), the sudden decay result differs 
from the non- instantaneous decay result by less than 1%. 
Hence, when constraining the curvaton model with the 
current observational constraints on /nl there is no need 
for an exact numerical calculation; using the sudden- 
decay approximation is sufficient. However, in the future 
experiments are expected to bring down the upper bound 
on |/nl|, and then constraining the curvaton model does 
require the numerical calculation presented here (or in 

Ha EH). 

In Fig. [3] we compare c/nl in the sudden-decay approx- 
imation with the non-instantaneous decay result. The 
sudden-decay result for <?nl is much more inaccurate than 
for /nl- However, the present observational constraints 
on (?nl are so weak that again the sudden-decay approx- 
imation may be sufficient. 

Let us end this subsection with a comment on numer- 



ical accuracy. Since the derivatives of N(p) in ([58)) 
(|60p involve subtraction of nearly equal numbers, the 
calculation must be carried out carefully The first re- 
quirement is that the integration step-size in N is small 
enough compared to Sp that the two initial values pi and 
Pi+i = Pi + Sp really lead to numerically different values 
for Ni = N(pi) and iVj+i = N(pi+i). The accuracy of 
N must be good enough to maintain enough significant 
figures in SN = JVj+i — Ni. The smaller the steps in Sp 
that we want to take the higher the accuracy in N that 
we need. We calculate the first derivative at p = p, as an 
average of two nearby gradients 



dN 
dp 



p=Pi 



1 / Nj-Nj-! 

2 V Pi - Pi-i 



N z+1 - N 
Pi+i - Pi 



(71) 



We use the same algorithm for the second derivative 
(with N replaced by the result of the calculation of 
dN/dp) and for the third derivative (with TV replaced 
by the result of the calculation of d 2 N /dp 2 ). As a result, 
the first derivative picks up contributions from 3 nearby 
points, the second derivative picks up weighted contribu- 
tions from 5 nearby points and the third derivative from 7 
nearby points. This procedure smooths out any residual 
numerical noise. 



C. An analytic approximation to the numerical 
result 



The analytic approximation of r(p), Eq. (|50| . with help 
of ([57]) gives 



2 x 1.24(1 -r) 1 - (1 -r) 



(72) 



Hence, from (|64[) we find an analytic approximation to 
the non-linearity parameter 



9j£_ 
9' 2 



(73) 



5W 
4r 2 I 



-2r + 2 x 1.24(1 - 



(1 - r) 1 - (1 - r)T^4 | 



The difference from the numerical result is non-negligible 
only when /nl is extremely close to zero. Indeed, we find 



NLfit - /NL 



/nl 



< 1% 



(74) 



if r < 0.501 or r > 0.542. The difference is larger than 
5% only when r G [0.528, 0.534]. 



V. PROBABILITY DENSITY FUNCTION 

Thus far we have calculated the second- and third- 
order corrections to the curvature perturbation produced 
by the curvaton decay from which the leading order terms 
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to the bispectrum and trispectrum can be calculated. 
However the 5iV-formalism allows us to describe the full 
non-linear probability density function on large scales for 
the non-linear primordial curvature perturbation defined 
in Eq. ©. 

Assume we have two random variables y and z, and 
the pdf of y is f(y). Furthermore, assume that the func- 
tional dependence of z on y is known, z = z(y), and this 
mapping is a bijection. Then the probability of z being 
in the interval z 2 ) is given by 



P(zi < z < z 2 )= I 

JZl 



rZ2 


dy 




dz 



f(y)dz. 



(75) 



where the absolute value is needed in the case that y(z) 
happens to be a decreasing function. Hence the pdf of 
the random variable z is 



/(*) 



dy{z) 



dz 



f[y(z)] 



(76) 



In the multi-variable case the derivative would be re- 
placed by the Jacobian determinant. For a Gaussian 
random variable, y, with mean fj, y and variance cr 2 , we 
have 



f(y) 



i 



27rcr2 



(77) 



Since the first-order primordial curvature perturba- 
tion, Ci, depends linearly on the initial Gaussian field 
perturbation, 8\*i we can take (i as our Gaussian "ref- 
erence" variable with mean fj,^ = and variance 



2 9' 

3 g Sx ' 



(78) 



With this goal in our mind we have already written 
all our non- linear expressions for £ as a function of £i- 
In the sudden-decay approximation we found an ana- 
lytic functional dependence C = C(Ci) and in the non- 
instantaneous decay case this function was found numer- 
ically. The mapping from £i to C is n ot actually a bi- 
jection as there can be several values of Ci which are 
mapped onto the same value of (. See Appendix A for 
the sudden-decay case. Calling these values £i •, it is now 
easy to calculate the pdf of the non-linear primordial cur- 
vature perturbation 



/(c) = E 



f g (Cu), 



(79) 



Ci=Ci 



where f g (Ci) it the Gaussian pdf with fi = and a = uq 1 . 

For simplicity we will assume in the rest of this section 
that there is no non-linear evolution of the curvaton field 
before it begins to oscillate, so that g cx x*, i- e -j 9^ n ' = 
for n > 2. In principle one could also carry through 
the non-linear evolution of the curvaton field into the 
full numerical calculation of the pdf for the primordial 
curvature perturbation. 



At the end of Appendix A we derive /(£) in the sudden- 
decay approximation, see Eq. (|A14[) . Here we continue 
by demonstrating the calculation of /(£) up to second 
order, which we will call /2(C)- From (HJ) we have £ = 
Ci + |/nlCi u P to second order, i.e., 



-1 ± y/T+ 12/ N lC/5 



6/nl 

Substituting this into (|79|) we find 



(80) 



/2(C) 



1 



1 



v/1 + 12/ NL C/5 

x ^ e 4^(-i±\^W)] 2 /K) i (81) 
± 

if C > ~5/(12/nl), and /2(C) = otherwise. If we want 
to evaluate ji (C) for a particular value of the parameter r 
we can substitute /nl from (|41[) into (JHTJ) in the sudden- 
decay case, or use /nl( - '') as shown in Fig. [5] in the non- 
instantaneous decay case. 

In Fig. IDJa) we compare the fully non-linear pdf /(£) 
(non-instantaneous decay) and /sd(C) (sudden decay) to 
the Gaussian / fl (Ci) in the case when there is large non- 
Gaussianity (r = 0.00028, p = 0.00030, / NL = 4432). 
In Fig. HJa) we plot the pdfs in the case when /nl has 
exactly the WMAP3 upper limit value (r = 0.010758, 
p = 0.011560, /nl = 114). When the non-linearity pa- 
rameter is very large, this kind of visual comparison re- 
veals the non-Gaussianity, but already with /nl = 114 
the / is virtually indistinguishable from the Gaussian f g . 
However, in Fig. BJb) and Fig. EJb) we plot f/f g which 
reveals the non-Gaussianity even when /nl = 114. 

A. Moments of the distribution 

The non-Gaussianity can be described quantitatively 
by calculating the moments of the pdf. Any pdf f(z) 
should give a unit total probability 

f(z)dz = l. (82) 
The mean can be calculated as 

Mz = / zf(z)dz . (83) 
and the i th moment m z (i) is defined as 

m,(i) = / (z-fiyf(z)dz. (84) 



The second moment is the variance {a 2 z ), the third mo- 
ment is called skewness, and the fourth moment kurtosis. 
As these moments can be extracted from the CMB maps 
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Figure 4: (a) Pdfs at r = 0.00028 (p = 0.00030, /nl = 4432). 
Red dashed line is for the sudden decay, /sd(C)> an< i black 
solid line for the non-instantaneous decay, /(C)- The solid 
green/grey line is the Gaussian "reference", / s (Ci)- (b) The 
ratio of non-Gaussian pdfs to the Gaussian one. 



Figure 5: (a) Pdfs as in Fig. H but now at r = 0.010758 (p = 
0.011560, /nl = 114). In this figure the Gaussian reference 
(solid green/grey line) is completely indistinguishable from 
the fully non-linear ("non-Gaussian") pdfs. (b) The ratio of 
non-Gaussian pdfs to the Gaussian one. 



it is enlightening to calculate the curvaton-model predic- 
tion for them. 

For Gaussian pdfs any odd moment (with i > 3) is 
zero, since the probability density is symmetric around 
the mean. The even moments of a Gaussian distribu- 
tion are easy to calculate employing partial integration 
to give m(4) = 3cr 4 , m(6) = 15cr 6 , m(8) = 105a 8 , 
m(10) = 945a 10 , m(T2) = 10395er 12 , etc. Any departure 
from these values indicates that the pdf is non-Gaussian. 
If odd moments differ from zero, there is an asymmetric 
deviation from Gaussianity. If even moments are smaller 
than in the Gaussian case, the pdf is more sharply peaked 
than the Gaussian. If even moments are larger, the pdf 
is wider. The set of moments {/j,,m(i)\i = 2.. .00} en- 
codes the same information of non-Gaussianity as our 
fully non-linear C(Ci) ( or the expansion £ = X^^Li Cn/n\). 
However, it should be noted that giving the value, for ex- 
ample, for TOf(3) is not simply equivalent to giving the 
value for /nl, because the moment picks up contributions 
from the fully non-linear £, not just from (1 + |/nlCi ■ 

It turns out that the moments can be calculated very 
accurately using the ^iV-formalism, even in the non- 
instantaneous decay case, since we need not to calculate 



the derivatives of the local expansion, N, as was done in 
calculating /nl or (?nl- At first it seems that we would 
need /(C), which includes a numerical derivative <i£i/d£: 

m c (i) = |(C-M)7(CR, (85) 

but substituting / from Eq. |79|) we end up with 

"*(*) = £ /(C-M)7 9 (Cy)dCy- (86) 

3 

Here the only numerically calculated quantity is £ (and 

Calculating the moments in the second order expan- 
sion and comparing to the results of fully non-linear 
calculation ([56]) we can address the question whether 
J2^=3 Cn/n\ gives an important contribution to £ and 
hence to the non-Gaussianity or whether the second or- 
der expansion £ « C«/ n ' is indeed accurate enough. 
To this end let us calculate in the second-order expansion 



11 




Figure 6: Third, fourth, fifth and sixth (from top to bottom) moments of the probability density function of the primordial 
curvature perturbation £ as a function of the linear transfer parameter, r. The predictions using the fully non-linear sudden- 
decay approximation are shown by the dashed red line (dot-dashed for negative values), and the second-order results are shown 
by the solid blue line (dotted for negative values). The fully non- linear numerical results are indicated by black crosses (stars 
for negative values). 



12 



the mean 



/ C/ 9 (CiRi - / ( Ci + |/nlCi 2 ) / s (GRi 



g/NLCT Cl 



(87) 



the variance 



(^1 + jj/NLCl'-^/NL^ /„(Cl)dCl 



^+2-/^4, 



the skewness 



ma(3) = 6 (|/ NL ) +8 Q/ NL ) 



(88) 



(89) 



and the kurtosis 



m 2 (4) = 34 +60 ^+60 (|/nl) <■ (90) 

For higher moments in the second-order expansion we 
find 

m 2 (5) = 600/^^+680 0/^ a\ 



544 -/nl 4 



(91) 



2 x 10 5 steps integrating TV from the Friedmann equation, 
and we calculate £ for 6001 equally spaced values of d 
in the range [—6 x 10~ 3 , 6 x 10~ 3 ]. Thus to produce a 
pdf for a fixed value of r we need 10 9 integration steps.) 
We compare the fully non-linear non-instantaneous decay 
result to the second-order results and to the fully non- 
linear sudden-decay result. 

We find that results obtained using the fully non- 
linear sudden-decay approximation agree well with those 
obtained from the full non-instantaneous decay. The 
sudden-decay approximation accurately predicts the mo- 
ments of the distribution for small values of r where the 
non-Gaussianity is largest, and only fails to give the pre- 
cise values of r where the moments cross zero, very simi- 
lar to what was seen previously when evaluating the non- 
linear parameters /nl and #nl m Sec. IV. 

The expressions for the moments, given in Eqs. (|87f - 
22), calculated using only terms up to second-order in 
perturbation theory [but using the full numerical value 
f° r /nlM] are an excellent description of the odd mo- 
ments of the distribution for all r. For even moments, 
the second-order expressions give the correct order mag- 
nitude, but cannot always reproduce the correct sign of 
the moments for r > 0.1. In particular we see that the 
even moments of the distribution predicted at second- 
order are always larger than the Gaussian value (setting 
/nl = 0), whereas the full numerical results show that 
the even moments can be less than the Gaussian value. 
To describe the deviations from Gaussianity in the even 
moments we need to include third-order terms. For ex- 
ample the variance, to third order, is given by 



and 



m 2 (6) = 154 + 1170 ( -/ NL ) < 



+ 9060 0/nl) 4° + 6040 0/nl) <jf t ■ (92) 

From Eqs. (f8"9")) and (|9"Tj) we find the leading order pre- 
diction 



m 2 (5)/<7 5 



= 10, 



m 2 (3)/a3 

and from Eqs. (|9T))) and (f9"2"| we expect 
m 2 (6)/cr 6 - 15 1170 



m 2 (4)/(T 4 - 3 



60 



19.5, 



(93) 



(94) 



at least when /nl is small. 

In Fig. [6] we plot the moments from the third up to 
the sixth one as a function of r. (For each value of r, 
it takes couple of hours to calculate the moments in the 
non-instantaneous decay case with our code on a typical 
PC. This comes about because we want to find N with 
a relative error of less than 10 -20 in order to have a 
sufficiently accurate result near to the peak of the pdf 
where C = 5N is extremely close to zero. We need about 



2(1/ 



15 ( 



2^5nl) 



NL 



6 (4<?nl) 



(95) 



where the expression in square brackets gives the leading 
order correction to the Gaussian result. This correction 
can be negative due to negative <7nl when /nl is small. 
On the contrary, the skewness up to third order is 

m 3 (3) = 6 (|/ NL ) < + [8 (|/ NL ) 3 + 72 (|/ NL ) 

i + 270 (f/ NL )(igNL) 2 <, (96) 



i2 %5nl) 



where the first term is the leading order correction to 
the Gaussian result [m(3) = 0]. As seen this correction 
does not have i?nl term so that already the second order 
expansion leads to approximately correct results. 

Although not shown in the Fig. [6] we have verified that 
the third-order perturbation theory (using the numerical 
results for /nl and (?nl) accurately reproduce all the mo- 
ments of the distribution as a function of r at least up to 
and including the sixth moment. 

The results (see Fig. [6]) obey the predictions of 
Eqs. (|93|) and ([94]) . In particular, the fully non-linear 
numerical calculation also reproduces the predicted ra- 
tios of the moments. 
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VI. VARIANCE ON SMALL SCALES 



and thus wc have 



We now consider the effect of a (possibly) large contri- 
bution to the curvaton density from smaller scale modes, 
compared with the cosmological scales probed directly, 
for instance, by the CMB anisotropics. This situation 
was rece ntly discussed by Linde and Mukhanov [3l[ (see 
also [13, l43| ) . Such smaller scale modes might contribute 
significantly to the average curvaton energy density on 
larger scales if the curvaton field power spectrum rises 
on smaller scales, or if some of the fraction (even if it is 
a small fraction) of the energy from the inflaton decay at 
the end of inflation is transfered to the curvaton [31( . In 
cither case wc will describe this by a small scale variance 
in the curvaton field up to some averaging scale 



A 



(Sx 2 



(97) 



The key observation is that these small-scale field fluc- 
tuations on spatially-flat hypcrsurfaces are uncorrelated 
with the field perturbations on larger scales. Thus there 
is an additional contribution to the average curvaton en- 
ergy density 



1 



Px 



m 2 (l + A 2 s ) x 2 



(98) 



which is homogeneous on large scales. In effect the cur- 
vaton density can be spilt into two parts: one that is 
perturbed on large scales, and one that is not. 

We can include the contribution from this small-scale 
variance in our non-linear expression (|34j) for the curva- 
ture perturbation, £, in the sudden-decay approximation 
where the curvaton decays on a uniform-density hyper- 
surface, to give the equation 



(1 - £! Xi doc)e 



-K 



1 



1 + A 2 x 
A? 



1 + A 2 x 



^ x ,dc C e- 3C = 1 , (99) 



where we have set to zero any pre-existing perturbation 
in the radiation, £ r = 0. 

At first order this shows how the resulting curvature 
perturbation on large scales is suppressed by the small- 
scale variance: 



Ci 



1 + A : 



2^X1 



(100) 



where r is given by Eq. (|38|) in the sudden-decay approx- 
imation, and we use £ x i to denote the fractional field 
perturbation on large scales, given in Eq. (j2"3")h 

However small-scale variance also affects the non- 
Gaussianity at second- and higher-orders. At second- 
order we find 



3(1 + A?) 



2r 



99 



-2-r 



(101) 



/NL = (1 + Aj)-(1 



99 

g' 2 



5r 
~6~ 



(102) 



Note that the small-scale variance, although homoge- 
neous, is not equivalent to additional homogeneous ra- 
diation due to its non-relativistic equation of state. 

If we allow for non-instantaneous curvaton decay we 
find 



4c 



5 d 2 N/dp 2 

6 (dN/dp) 



2 ' 



(103) 



where any A 2 dependence in the second term on the 
right-hand-side cancels out so that we can use the numer- 
ical results presented in Sec. IIVI to evaluate this term. 



A. Observational constraints on small-scale 
variance 

The fact that the non-linearity parameter grows with 
the small-scale variance means that we can constrain the 
small-scale variance from constraints on the non-linearity 
parameter /nl on larger scales. 

In practice the non-linearity at each successive order 
still depends on the non-linear evolution function for the 
curvaton field, g(x*)- Hence we cannot rule out mod- 
els where the small-scale variance is large, but its effect 
is precisely cancelled by the non-linear evolution. For 
simplicity we assume in the following that the non-linear 
evolution is negligible so that g" and higher derivatives 
can be set to zero. 

Recalling, that r < 1 and — 54 < /nl < 114 (from 
WMAP3 0]) we find an upper bound 



At < 90 . 



(104) 



B. 



Observational constraints on variance on CMB 
scales 



On the scales directly probed by CMB observations, 
the constraint on the variance will be much tighter, since 
in addition to /nl constraint we observe 



(Cl)cMB = A 2 , 



(105) 



with A 2 ~ 6.25 x lO" 10 . Substituting (fT00|) with j23]) 
into the left hand side we get 



4 A 2 

4 CMB 2 _ a2 
9 (1 + A 2 ) 2 



where 



A 2 

"CMB 



(*x 2 : 



CMB 



(106) 



(107) 
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Equation (I106[) gives 
A 2 



9, 2 (1 + A 2 ) 2 

CMB ~ ~ A A ~2 ' 



(108) 



and eliminating (1 + A 2 ) 2 /r 2 with help of (|103|) we end 
up with 



A 



CMB 



9 2 / 4 2 d 2 N/dp 2 
-,A -/nl - - i 



4 15' 3(dN/d P y 



(109) 



The maximum of the absolute value of the second term 
in the parenthesis is numerically found to be always less 
than 2. Thus, employing the triangle inequality, we find 



A C mb<^ 2 (|^nl|+2 



(110) 



But the WMAP3 upper limit for ||/nl| is 91, which im- 
plies 



A CMB < x 932 = I- 2 x 10_5 • ( m ) 



VII. CONCLUSIONS 

In this paper we have presented for the first time the 
fully non-linear probability density function (pdf ) for the 
primordial curvature perturbation on large scales in the 
curvaton scenario using the 5A-formalism. By solving 
the non-linear evolution equations in an unperturbed 
(FRW) universe one can construct the local expansion 
upto a final uniform density as a function of the initial 
curvaton field value, N(x)- Assuming a Gaussian form 
for the initial field distribution on large scales (as would 
be expected for a weakly coupled scalar field after in- 
flation) it is straightforward to construct the probability 
density function for SN and hence the non-linear curva- 
ture perturbation £, defined in Eq. @. This procedure is 
particularly simple in the case where the local expansion 
is a function of a single scalar field, such as the curvaton, 
but it is also straightforward to apply to multiple fields 
whose initial distributions on large scales are known. 

In the sudden-decay approximation where it is as- 
sumed that the curvaton decays instantaneously, when 
H ~ r, we have presented a simple non-linear analytic 
expression, Eq. (f34|) , relating the primordial curvature 
perturbation to the initial curvaton perturbation. We 
have compared analytic results in the sudden-decay ap- 
proximation, with our results derived from direct numer- 
ical integration of the full coupled equations for the local 
radiation and curvaton energy densities and found good 
quantitative agreement. 

In particular we have calculated the leading-order con- 
tributions to the primordial bispectrum and trispectrum, 
including for the first time the effect of third-order terms 
in the curvature perturbation. In some cases [35l.l3l| non- 
linear evolution of the curvaton field on super-Hubblc 



scales, after Hubblc-cxit during inflation, but before the 
curvaton begins to oscillate about the minimum of its 
potential, could lead to a suppression of the leading or- 
der contribution to the primordial bispectrum. We have 
shown that in this case there will instead be a large con- 
tribution to the primordial trispectrum, unless there is 
an additional cancellation in the third-order term. 

We have computed numerically the moments of the 
pdf for the primordial curvature perturbation up to and 
including the sixth-order moment for a range of values 
of the linear transfer coefficient, r. To accurately repro- 
duce the even moments of the distribution we need to 
go beyond the second-order terms in the curvature per- 
turbation (described by the non-linearity parameter /nl) 
and include higher-order terms. 

One example of how non-Gaussianity can be used to 
constrain model parameters is the case when the cur- 
vaton field has a large variance on small scales, as re- 
cently proposed by Linde and Mukhanov [3l|. In this 
case the suppression of the linear curvature perturbation 
is accompanied by an increase in non-Gaussianity. We 
have shown that in this case limits on the primordial bis- 
pectrum can be used to place limits on the small scale 
variance. 

The calculations presented here should enable the cur- 
vaton model to be subjected to a range of tests of non- 
Gaussianity, going beyond just the bispectrum. In the 
simplest models (neglecting non-linear evolution of the 
field before it decays) the non-Gaussianity is a function 
of a single parameter, r, which is the linear transfer co- 
efficient relating the first-order primordial curvature per- 
turbation with the curvaton perturbation at Hubble-exit 
during inflation. Multiple tests of the form of any pri- 
mordial non-Gaussianity could offer consistency tests of 
the curvaton scenario. 
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Appendix A 

In this appendix we solve the primordial curvature per- 
turbation ( as a function of initial Gaussian field per- 
turbation Six in the sudden-decay approximation. We 
also solve the inverse problem, i.e., find 5\x/x ( or Ci) as 
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a function of £. Using these results we derive an ana- 
lytic expression for the (non-Gaussian) probability den- 
sity function of £; /(C)- 

We can rewrite Eq. |35|) in the form 



4r 



3r-3 



= 0, 



(Al) 



where r = 3f2 X) dec/(4 — f^ x ,dcc)- This is a fourth degree 
equation for X = e^. The solution of this full non-linear 
equation which gives the primordial curvature perturba- 
tion as a function of initial Gaussian curvaton field x* ( x ) 
is 



Hence (|A7|) simplifies to a second degree equation for 
5\x/X- The solutions are 



X 



= -1 ± 



3 + r 
4r ' 



,3C 



3r- 3 
4r 



-C 



-,1/2 



(A9) 



where the "+" sign corresponds to a small perturbation 
and "— " sign would give |<$ix/x| ~ 1- A n alternative 
Gaussian "reference variable" is the linear end result Ci . 
From J57J) and (J23J) we have Ci = In Sec. [V] we 

will need a derivative of this Gaussian random variable 
Ci with respect to C- Using (|A9[) we easily find 



with 



where 



X = K 1 ' 2 



C = InpO , 
1 + y/ArK- 3 ' 2 - 1 



(3 



A d = f e 3 ^ = 



\l/3 



<?(x*«) 



ff(x* 



A" 



def 1 



pl/3 + ( r _ l)( 3 + r )l/3p-l/3] 

P = (Ar) 2 + [(Ar) 4 - (3 + r)(r - 1) 



31 1/2 



(A2) 

(A3) 

(A4) 

(A5) 
(A6) 



The inverse problem (solving initial xlx as a function 
of C) is much simpler. Namely, Eq. (|A1[) gives immedi- 
ately 



4r 



-e 3 < + 



3r-3 
< 

4r 



(A7) 



and here e 3C * = 5 2 [x*(x)]/5 2 (x*) = % 2 (x)/x 2 . 

Assuming that there is no non-linear evolution be- 
tween the Hubble exit and start of curvaton oscillation 
Lg (rl) = for n > 1] the left hand side of (|A7|) is exactly 
[see Eq. (fl9j) ] 



1 + ^) =l + 2^+(^) • (A8) 
I 



<*Ci± 



1 



4r 
3 + r 



3C 



.3C . 



3r-3 
4r 

3r- 3 



-C 



-1/2 



. (A10) 

4r 4r J 

Hence the full non-Gaussian probability density function 
for C is 



/(C) = /-(C) + /+(C), 



where 



and 



/±(C) 



/ s (Ci±) 



/ 9 (Ci±), 



-C?±/(2^) 



2^ 



(All) 



(A12) 



(A13) 



with Ci± being 2r/3 times the rhs of Eq. (|A9|) . Here 
fg(Ci) i s the Gaussian pdf for the first order perturbation 
Ci with variance c 2 = cr 2 ^ (~ 6.25 x 10~ 10 to match the 

observations), and mean /i = 0. In practice, /-(C) could 
be neglected, because / s (Ci_) is typically of the order 
cxp(— 10 10 ). Substituting all ingredients into (jAlip the 
pdf reads 



/sd(C) 



l l r 



2ttct2 2 



(3 + r)e 3C + (1 -r)e~ c 



3 + r o,* 3r — 3 r 
e 3C H e~ c 



-1/2 



-1 ± 



3 + r 3r — 3 



4r 



4r 



1/2' 



'(2a 



(A14) 



where the subscript SD reminds us that this is the sudden-decay result. 
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